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ABSTRACT 

Context. Brown and Mallik (BM) recently showed that, for hot sources, recombination of non-thermal electrons (NTR) onto highly 
ionised heavy ions is not negligible compared to non-thermal bremsstrahlung (NTB) as a source of flare hard X-rays (HXRs) and so 
should be included in modelling non-thermal HXR flare emission. They further claimed that, in some cases, NTR can be much larger 
than NTB with important consequences for flare physics. 

Aims. In view of major discrepancies between BM results for the thermal continua and those of the Chianti code and of RHESSl solar 
data, we critically re-examine and correct the BM analysis and modify the conclusions concerning the importance of NTR. 
Methods. The BM and Chianti element abundances and ionisation fractions as a function of temperature T for relevant elements are 
found to agree well. Although the analytic Kramers expression used by BM is correct for the purely hydrogenic recombination cross 
section, the heuristic expressions used by BM to extend the Kramers expression beyond the 'bare nucleus' case to which it applies 
had serious errors. BM results have therefore been recalculated using corrected expressions, which have been validated against the 
results of detailed calculations. 

Results. The BM results are found to be correct for NTR onto Fe 26+ and a factor of 2 too high for Fe 25+. Thus, at high enough 
T > 40MK for these to exist, such NTR strongly dominates NTB in the deka-keV range just as BM claimed. However, at such T, 
thermal continuum dominates NTR and NTB in this energy range unless the non-thermal electron density is a very large fraction, f^, 
of the total as in some coronal HXR sources. At T a 10-30 MK the dominant ions are Fe 22+, 23+, 24+ for which BM erroneously 
overestimated NTR emission by around an order of magnitude. Contrary to the BM claim, NTR in hot flare plasmas does not dominate 
over NTB, although in some cases it can be comparable and so still very important in inversions of photon spectra to derive electron 
spectra, especially as NTR includes sharp edge features. 

Conclusions. The BM claim of dominance of NTR over NTB in deka-keV emission is incorrect due to a serious error in their analysis. 
However, the NTR contribution can still be large enough to demand inclusion in spectral fitting, the spectral edges having potentially 
serious effects on inversion of HXR spectra to infer fast electron spectra. 
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1. Introduction 

iBrown & Mallik' ('2008', '2009^ (BM) emphasised the importance 
of non-thermal recombination (NTR), previously neglected as a 
mechanism for production of solar hard X-rays (HXRs). BM ob- 
tained analytic Kramers approximations for recombination and 
bremsstrahlung continua and showed the resulting spectra from 
Maxwellians and cut-off power-laws to elucidate the importance 
of NTR. BM also used typical flare parameters to support their 
findings and concluded that NTR is not only an important com- 
ponent of the HXR spectrum in solar flares, but can be the dom- 
inant source in hot plasmas. This had serious implications for 
spectral inference of electron distributions and for electron num- 
ber and energy budgets. However, BM did not compare their 
results w ith actual flare HXR data from instruments such as 
RHESSl (Li n et aljr200 2') nor their thermal spectra wi th com- 
puted atomic/radiation databases such as Chianti (e.g. De re et all 
l20~09l) . In this paper, we discuss such comparisons, which re- 
vealed serious discrepancies overlooked by us till now and which 
led us to discover blunders in BM that we amend here, with 
apologies. 
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Chianti includes only thermal emission (bremsstrahlung 
(TB) and recombination (TR) plus lines) and not non-thermal 
bremsstrahlung (NTB) nor NTR, so we can only compare it with 
BM results for the Maxwellian case. Fig. [T] re veals that the 'ex- 
act' Chianti and the approximate BM results agree well for TB 
spectra but that the TR emission BM predicts is much higher 
around the 10 keV range due to the Fe and Ni emission in BM 
which shows large r ecombination edges (around 7-8 keV) almost 
invisible in Chianti. ILandil (l2007h emphasised the importance of 
TR compared to TB, but his predictions and spectra were also 
quite different from BM. (Unlike Chianti, BM spectra did not 
include emission lines but even the lines in Chianti are smaller 
than BM's recombination edges). After BM, we also recognised 
that Chianti thermal spectra can fit pretty well to what RHESSl 
observes, a fact pointing strongly toward some flaw in the BM 
modelling. These discrepancies motivated this paper which un- 
covers and corrects several errors in BM. 

Having checked all our algebra and arithmetic, co mpared 
our Kramers purely hydrogenic expressions with those of lLandil 
(1200 7) and others and checked the accuracy of Kramers' cross- 
sections against more exact atomic results and also of BM 
against more exact Maxwellian-weighted emission rates (all 
good within about 15-20% at worst), and that BM element abun- 
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dances and ionisation fractions were in hne with Chianti, we had 
to ask whether there was a more basic flaw in the BM heuristic 
extension of Kramers expressions for H-like ions to He- and Li- 
like and Li-like to Be-like. We found that there was. A crucial 
clue lay in a closer look at Fig.[T] which shows that the energies 
and sizes of edges for all the ions in the BM code agree quite 
well with those in Chianti, except for Fe and Ni, which proved 
discrepant. These are the only significant contributors which are 
not near full ionis ation. 

Sect. 2.3 of iBrown & MaUikI (|2008|) stated "Noting that 
Qr oc 1 /m^ we include here only recombination to m = 1 (in 
the sense of the lowest empty level of the ion - hydrogenic with 
Z = ZefF - not of the atom). Higher m contributions are weaker, 
being oc 1/m^ though extending to lower energies with edges at 
■^efF^/'"^ "' where is the recombination cross-section. Herein 
lies the source of the BM blunder and the source of discrepancy. 
The relevant 'm' in the cross section for the lowest unfilled level 
is not m = 1 in general but the principal quantum number for 
that lowest unfilled level of the ion involved. For Fe26-i- and 25+ 
m - \, but for Fe 24+ to Fe 11 +, for example, m + I but - 2! 
So the BM interpretation of m was at fault. Henceforth, we shall 
denote the principal quantum number by its standard notation 

to avoid further confusion. In ionic species which already 
have two or more electrons present, the smallest n value nmin is 
2, and if there are 10 or more electrons, it is Wmin - 3 and so 
on. (Recombination rates to levels with n > Wmin still fall off as 
so are rather small in comparison). The consequences of 
this BM error are : 

- since, for typical hot flare temperatures of 20-30 MK, Fe 24+ 
and Fe 23+ are the most abundant Fe ions, n^i„ = 2, not L 
Consequently, the magnitude of the recombination emission 
is down compared to BM by a factor around l/n^^^ = 1/8. 

- equally importantly, the locations of recombination edges 
for these Fe species are no longer at Z'^^x^ or ^lf(X + 

the presence of a low energy cut-oflF E^, but at Z^^x/^ or 
Z^gX/^ + Ec respectively. 

- a separate issue missed by BM is that the Kramers formula 
assumes, and is applicable to, recombination into an empty 
shell. For highly-charged ions, reasonable account is taken of 
initially partially filled n-shells by applying a 'vacancy fac- 
tor' p„ to the usual Kramers formula. The simplest choice 
for p„ is Ny/N,, where N„ = 2n^ is the total electron occu- 
pation number of an n-shell and A^,, the number unoccupied. 
For recombination of a H-like ion to « = 1 of a He-like ion, 
pi - 1/2. For recombination into a partially filled n - 2 
shell, A^2 = 8 would follow. However, a more accurate re- 
sult is obtained on recognising that recombination into n = 2 
in the Kramers formula is dominated by the 6 p-states, i.e., 
Uttle of the Kramers n - 2 result arises from recombina- 
tion into s-states, at least at the electron energies of interest 
here. Thus, we take N2 = 6 and A^,, the number of unoc- 
cupied 2p states, i.e. pi - ^ for Li and Be-like initial ions 
and p2 - 5/6, 4/6, ... 1/6 for B- through F-like initial ions. 
Comparisons of such modified Kramers cross sections have 
been made with the results of detailed calculations using the 
AUTOSTRUCTURE code (cf. Badnell 2006) for initial H- 
like through to F-like Fe ions and agreement to within 20% 
is obtained at the electron energies of interest here. 

These considerations explain why all the elements apart 
from Fe and Ni produced edge positions and strengths in BM 
reasonably comparable with Chianti, since those elements are 
almost fully ionised and nmin = 1 at relevant temperatures. The 



high Zed ions Fe and Ni are diff'erent since they are predomi- 
nantly in ionic stages with n - \ filled and = 2. (For super 
hot plasmas, say > 40 MK, Fe 25+ and Fe 26+ with n - \ start to 
become important, but at those temperatures the thermal emis- 
sion begins to dominate emission up to 50-60 keV and to dom- 
inate any NTR emission and recombination edges in particular, 
an issue we discuss later). In this paper we address the implica- 
tions of these corrections. For full details see Mallik (2010). The 
BM results were mainly for the widely used, though unphysical, 
electron flux spectrum F{E) oc £ ?A E > Ec and zero below 
that. For generality of our corrections here we also give the NTR 
expression for general F{E). 

Sections 2 and 3 provide the key amended expressions and 
results to replace the erroneous ones in BM and new plots to 
emphasise the differences from BM, and to re-assess the impor- 
tance of NTR. Although the magnitude of our NTR predictions 
has been much reduced, it is still significant and requires consid- 
eration, especially in cases where the spectral index 6 is large, 
the low-energy cut-off Ec is small, and temperature T is between 
20 and 30 MK. Our overall conclusions are discussed in Section 
4. 



2. Alterations to relevant expressions 

Revisions of the recom bination cross-sec t ion ex pressions given 
inEqs. (11) and (12) of 'Brown & MaUikI (l2008l) have to be car- 
ried through to all other recombina tion expressio ns, e.g. Eqs. 
(16, 19, B.2, B.5, B.12) in IBrown & MaUikI (l2008h and Eq. (2) 
in Brown & Mallik (20()9h . to show the sum over n > Wmin, with 
«min 1 in general. So, th e revised thermal recom bination ex- 
pression (from Eq. (B.5) of iBrown & MaUikI (|2008|) ) as a func- 
tion of photon energy e with Vz^g — Z^^x is (all notation as in 
BM) 



XexpP^li- ; e>yz,,/«' 



■^R,hc,m(f) 



X ; e < Vz^Jn\ (1) 

The revised thin-targe t non-thermal recombin ation expression 
(replacing Eq. (B.2) of lBrown^ MaUikI (120081) 1 is: 



^R,„„(e) = (5-1) 



327r rlx^ 2npALFc v V >- 1 
2j Zj /^"tRZ.,r3 



V3a f El 



-S-l 

; £>Ec + Vz,Jn^ 
X ; e<Ec + Vz^Jn^ (2) 



for the cut-off power-law F(E). We also give the corrected thin 
target expression for general F(E) in Eq.[3]to enable readers to 
extend our results: 



327rr^ynpV 



Z «>nmin 



AzZ^ F(E-Z^X/n^) 
E-Z^Xln^ 



(3) 



The thick-target NT R expression (revised fromEq. (B.12) of 
IBrown & MalUkf j2008h ) is 



Brown, Mallik & Badnell: Non-thermal recombination revisited 



3 



Z^fr n>«min 



3V3a 



X ; e < yz,,/nl 



(4) 



Eqs. (|2]i and (|4|l also correct several typos in expressions 
(B.2) and (B.12) of BM respectively. Thick target results are 
for the total emission rates over continuously injected electron 
collisional lifetimes. 

Although all n values should strictly be included for all ions, 
in practice only the recombination to the lowest available level is 
significant, because the rates are oc Hence, in our calcula- 
tions, we use only n - - 1 for all elements up to Ca, while 
from Ca onwards, we have n = nmi,, = 2 for the appropriate ions, 
e.g. Fe 23+, Fe 24+ or Ni 25-I-. With very high temperatures Fe 
25-1- and 26-1- also become important. For these species, Wmin = 1 ■ 

The ionisation equilbrium used was a fit to the standard 
steady state coronal collisional ionisation, optical depths being 
negligible. In Table [T] we have listed the abundances and ^rz 
for Fe 25+ and Fe 26+ (Arnaud & Raymond 1992) for a set 
of four temperatures, and these can b e seen in additio n to th e 
relevant values in Tables 1 and 2 of iBrown & MallikI (l2008h . 
where these values are listed for fully ionised Fe (Fe 26+) at 
r » 100 MK and for the most abundant Fe species at 20 MK. 
Here we have also included the ratio of non-thermal free-bound 
(due to only Fe 25+ and 26+ respecti vely) to total free-free fl ux, 
R = Jntr^,26..25JJntb (fromEq. 2 of .Brown & Malhkl (I2009I) '). at 
e = 20 keV and (5 = 5 to illustrate at what temperatures Fe 25+ 
and Fe 26+ start becoming significant contributors to the NTR 
flux. 

Note that we use the variable fc to define the fraction of total 
electrons that are accelerated, i.e. /. = rtc/rip, where is the 
density of non-thermal electrons. 

3. Results 

Some of the plots from BM have been reproduced here to show 
the essential changes in our results. We have also added some 
new plots to show interesting new results. We have not repro- 
duced all the plots as that would be repetitive but, from these 
sample plots, it should be clear what the changes are and what 
they imply. Throughout this section we have used a resolution 
of 0. 1 keV bins, unlike in BM, where primarily 1 keV bins were 
used. We do this so that the finer details can be appreciated. 

Firstly, Fig. |2] shows our thermal spectra compared with 
those of Chianti at 20 MK. It is clear here that the major 
BM discrepancies from the plot in Fig. [T]have been removed. 
The slight dififerences between our model spectra and those of 
Chianti can be attributed to finer details such as the use of 
Kramers' cross-section versus more accurate calculated cross- 
sections. But while the thermal edges are much smaller, they 
can still form a significant proportion of the flux and may be 
detectable. 

In Fig. |3] we redo the bottom-left plot from Fig. 1 of 
iBrown & Mal Uk (2008), where ^ 10 keV, 5 = 5 and T = 20 
MK (hot), showing only the thin-target non-thermal spectra. The 




e (keV) 



Fig. 1. The original BM thermal model spectra (dot-dash blue is 
TB; solid red is TR) compared with Chianti's (dotted green for 
TB; dashed black for TR) for T = 20 MK. Courtesy E. Landi for 
the Chianti spectra. 
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Fig. 2. The revised thermal model spectra compared with the 
same Chianti spectra for 20 MK; the colour/linestyle-coding is 
the same as in Fig. [T] 



main point to note here is that at no point is NTR greater than 
NTB, but it can still be a significant contributor to the total non- 
thermal flux in and around the cut-off energy, in this case 10 keV. 
However, the thermal emission at these energies is quite likely 
to limit the detectability of the non-thermal edges for these tem- 
peratures. Nevertheless, it is important to note that all the major 
edges occur at around 2-4 keV (or 2-4 keV above the cut-off in 
the case of NTR thin-target) because our revised Vzirp- value for 
the important species of Fe, namely Fe 23+ and Fe 24+ in this 
case, is about 2-3 keV. So these important Fe edges coexist with 
edges from other elements like Si and O that have significant 
values of (rz and whose Vz - 3-4 keV. These are almost fully 
ionised, so nmin = 1 for them. 

Fig.|4]shows the non-thermal thin target NTR:NTB ratio for 
the same parameters and can be contrasted with the bottom-left 
plot of Fig. 2 of Brown & MalHk (2008). These plots are vir- 
tually identical in shapes, though the revised plot has smaller 
values for the ratio, as expected. 
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Table 1. Significance of Fe 25+ and 26+ at 4 different temperatures for e = 20 keV and 6 = 5. 
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Fig. 3. NTB (dashed blue) and NTR (solid red) spectra from the 
revised thin-target model with the solid black curve denoting the 
total non-thermal flux multiplied by 10. Ec = 10 keV, 6 = 5 and 
r = 20 MK. 




e (keV) 



Fig. 4. Thin-target NTR:NTB ratio for the same flare parameters 
as in Fig. [3] 



Fig. |5] is the corresponding modified thick target plot of 
Figure 5 from Brown & Mallik (2008), and is our predicted 
spectrum for the 2002 April 14 flare parameters. Here too the 
non-thermal edge is comparable to the total non-thermal flux at 
the cut-off energy of 20 keV and might be detectable in observed 
spectra. The fact that NTR can be comparable to NTB for these 
particular real flare parameters is of particular significance. Also, 
the sharp nature of the observed edges can have a large effect on 
spectral inversions, an issue we are currently investigating. 
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Fig. 5. Our revised model thick-target plot spectra for the 2002 
April 14 event. The dashed blue curve is TB, dotted red is TR, 
dot-dash green is NTB and solid black is NTR. The solid ma- 
genta curve is the total flux multiplied by 10; small TR and NTR 
edges are clearly vis ible but much less distin ct from what was 
predicted in Fig. 5 of lBrown~& Mallikl(l2008h 



Since Fe 25+ and Fe 26+ start becoming important only for 
temperatures above 30 MK (Table [1]), we decided to look at our 
model spectra for each of the four different temperatures while 
keeping = 10 keV and 6 = 5 fixed. With fc = 0.1, we plot re- 
sults in Fig.|6]for the thin-target case. Although the Fe 25+/26+ 
edge becomes more prominent for the higher temperatures and 
NTR can be greater than NTB at these edges, higher tempera- 
tures also mean an increased thermal emission (where TR be- 
comes less and less important than TB), which tends to hide 
these non-thermal features. Hence, if NTR is to be observed, 
the optimum temperature would be about 20-30 MK. The same 
temperature range is also optimum to observe TR edges, as is 
validated by Landi (2007), who also predicts that TR starts be- 
coming less significant for T > 25 MK. 

The surface plot in Fig. Q is analogous to Fig. 2 of 
iBrown & Mallik (2009), showing the regions of relevant impor- 
tance. It is evident that NTR is never dominant but still can be 
several 10s of percent of the NTB flux. Once again, increas- 
ing fc to 0. 1 pushes the thermal-dominated line by about 5 MK 
to the right, but of course never changes the NTR: NTB ratio. 
However, in the (e, T) range of 20-30 keV and 20-30 MK re- 
spectively, NTR can still be significant enough for its inclu- 
sion in spectral analysis packages to be essential, moreso when 
fc is even larger as that increases the region of importance. 
As an extreme case, we produce Fig. [8l where f. — > 1 (as 
in Kruckeretal. (2008)). Thus, when fc is large, NTR domi- 
nated regions still exist in (T, e), albeit at rather high T and over 
a modest e range. Nonetheless, R is substantial (> 0.5) for a 
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Fig. 6. The revised thin-target spectra for 4 different tempera- 
tures with fc = 0.1. Note that the non-thermal emission measure, 
EMc - fcEM, where EM - 2ALn^ is the thermal emission mea- 
sure. The colour/linestyle-coding is the same as in Fig.|5] 
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Fig. 7. Regime plot with our revised model showing the relevant 
areas of importance in the (e, T) domain for E^ - 10 keV, 6 = 5 
and/c =0.01 
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Fig. 8. Regime plot with our revised model showing the relevant 
areas of importance in the (e, T) domain for E^ = 10 keV, 6 = 5 
and /c = 1 . This shows that NTR may still be a dominant source 
of HXRs when approaches unity. 

large (e, T) regime of 15-40 MK and energy range 20-35 keV, 
which is crucial when con sidering albedo spectral corrections 
(lAlexander & Browiil |2002|) and also electron number and en- 
ergy budgets. Even if fc, T and 6 are not so high, NTR can still 
be up to a 20-40% effect. 

(Note that all the other parameters for all these plots have 
been maintained throughout, viz.: density tip = 1.5 x 10" cm"^, 
flare loop cross-section area A = 2 x 10'^ cm^ and loop half- 
length L = 2.5 X 10*^ cm, mirroring thos e observed for the 2002 
April 14 event dVeronig & Brownll2004 ). 

4. Conclusions and implications 

The upshot of our revised model is that although NTR is (con- 
trary to BM) seldom dominant over NTB and the total thermal 
emission, it can still be substantial enough for its inclusion to 
be essential. The recombination edges clearly seen in our model 
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spectra (Fig.|5]l, though small, could have a large impact on elec- 
tron spectra inferred by spectral inversion since this involves 
differentiation of the photon spectrum. In extreme parameter 
regimes, NTR can be a substantial or dominant source of HXRs 
and so a large consideration in electron number and energy bud- 
gets and hence acceleration efficiency problems. However, for 
the majority of hot HXR sources, NTR is unlikely to be domi- 
nant but still significant enough to demand consideration. 

Another enhancement to the NTR contribution is the possi- 
ble enrichment of Fe abundance in the corona compared to pho- 
tospheric levels by a factor even hig her than the factor of 3 we 
have used based on 'Feldman (1992). Factors of ~ 10 or more 
have been suggested by Feldman et al. (2004), which would in- 
crease our NTR predictions by another factor of 3. In that case 
even for 'normal' flare parameters, NTR would start pushing the 
NTR-dominant boundary. Clearly, much work needs to be done 
to see if the First lonisation Potential (FIP) effect really does en- 
rich low FIP elemen ts, such as Fe, to the levels mentioned by 
iFeldmanet all (12004 . 

Further improvements to our model can be made by using 
more accurate recombination cross-sections instead of Kramers. 
But our research has revealed that Kramers' cross-section is not 
a bad approximation and is accurate to within 20% of calculated 
cross-sections for highly-charged ions. Another refinement in 
our model would be to include recombination onto 2 or 3 higher 
levels apart from the lowest available level of the ion. This too 
would make our model more accurate by a few percent at the 
most. In any case, it will be necessary to include NTR in spec- 
tral analysis packages such as SolarSoft to get reliable inferences 
from flare data. 

Acknowledgements. We gratefully acknowledge the financial support of a UK 
STFC Rolling Grant (JCB), of a Dorothy Hodgkin's Scholarship (PCVM), an 
ISSI Grant (JCB), and an STFC Grant (NRB). Discussions with E.P. Kontar, 
R.A. Schwartz, and E. Landi were invaluable in drawing attention to the fact that 
there was likely something amiss with the BM results, which we have addressed 
in this paper 

References 

Alexander, R. C. & Brown, J. C. 2002, Sol. Phys., 210, 407 

Amaud, M. & Raymond, J. 1992, ApJ, 398, 394 

Badnell, N. R. 2006, ApJS, 167, 334 

Brown, J. C. & Mallik, R C. V. 2008, A&A, 481, 507 

Brown, J. C. & Mallik, R C. V. 2009, ApJ, 697, L6 

Dere, K. R, Landi, E., Young, R R., et al. 2009. A&A, 498, 915 

Feldman, U. 1992, Phys. Scr, 46, 202 

Feldman, U., Dammasch, I., Landi, E., & Doschek, G. A. 2004, ApJ. 609. 439 
Krucker, S., Wuelser, J. -P., Vourlidas, A., et al. 2008, 12th European Solar 

Physics Meeting, Freiburg, Germany, held September, 8-12, 2008. Online at 

http://espm.kis.uni-freiburg.de/, p. 2. 84, 12, 2 
Landi, E. 2007, A&A, 476, 675 

Lin, R. R, Dennis, B. R., Hurford, G. J., et al. 2002, Sol. Phys., 210, 3 
Mallik. R C. V. 2010, PhD thesis. University of Glasgow 
Veronig, A. M. & Brown, J. C. 2004, ApJ. 603. LI 17 



